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We derive a semiclassical trace formula for the level density of the three-dimensional 
spheroidal cavity. To overcome the divergences and discontinuities occurring at bifurcation 
points and in the spherical limit, the trace integrals over the action-angle variables are per- 
formed using an improved stationary phase method. The resulting semiclassical level density 
oscillations and shell energies are in good agreement with quantum-mechanical results. We 
find that the births of three-dimensional orbits through the bifurcations of planar orbits 
in the equatorial plane lead to considerable enhancement of shell effect for superdeformed 
shapes. 

§1. Introduction 

The periodic orbit theory (POT)^^"^*^) is a nice tool for studying the corre- 
spondence between classical and quantum mechanics and, in particular, the inter- 
play of deterministic chaos and quantum-mechanical behavior. But also for sys- 
tems with integrable or mixed classical non-linear dynamics, the POT leads to a 
deeper understanding of the origin of shell structure in finite fermion systems such 
as nuclei, '^^'^^•'"^^^ metallic clusters, ^^•'"^^^ and mesoscopic semiconductors. "'^'^^"^^^ Bi- 
furcations of periodic orbits may play significant roles, e.g., in connection with the 
superdeformations of atomic nuclei,^^'^-*'^^-*'^^-*"^^-* and were recently shown to af- 
fect the quantum oscillations observed in the magneto-conductance of mesoscopic 
devices. ^^•''■^''^ This phenomenon is observed for some control parameters (shapes, 
magnetic field etc.) of the potential well, for which the orbits bifurcate and new 
type of periodic orbits emerge from the original ones. Examples can be found, e.g., 
in elliptic billiard and spheroidal cavity.^^'^^'^^^'^^^"^^^ In elliptic billiard, short dia- 
metric orbits with repetitions bifurcate at certain values of deformation parameter, 
and new orbits with hyperbolic caustics (butterfly-shaped orbit etc.) emerge from 
them. In spheroidal cavity, periodic orbits in the equatorial plane bifurcate, and new 
three-dimensional orbits emerge. 

The semiclassical trace formulae connect the quantum-mechanical density of 
states with a sum over the periodic orbits of the classical system. ^^"'^^ In these for- 
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mulae, divergences arise at critical points where bifurcations of periodic orbits occur 
or where symmetry breaking (or restoring) transitions take places. At these points, 
the standard stationary phase method (SSPM),*^ used in the semiclassical evalua- 
tion of the trace integrals, breaks down. Various ways of avoiding these divergences 
have been studied,^-*' ^^'^^^ some of them employing uniform approximations. ^^^^^^^ 
Here we employ an improved stationary-phase method (ISPM) for the evaluation of 
the trace integrals in the phase-space representation, which we have derived for the 
elliptic billiard^^^ and very recently for the spheroidal cavity.^^^ 

The singularities of the SSPM near the bifurcation points are due to the peculiar- 
ities of its asymptotic expansions. In the ISPM,^^)'^^) the catastrophe integrals^^^'^^^ 
are evaluated more exactly within the finite integration limits in the phase-space trace 
formula, and one can overcome the singularity problem due to bifur- 
cations, which occur when the stationary points lie near the ends of the integration 
region in the action-angle variables. We can also take into account the station- 
ary points outside the classically accessible region ("ghost orbits" ).^^ This method 
is particularly useful for integrable systems where integration limits are easily ob- 
tained. This theory has been developed in Ref. 28) for the case of the bifurcations 
through which periodic orbit families with maximal degeneracy emerge from the or- 
bits with smaller degeneracy. The essential difference between our method presented 
in this paper and that with the uniform approximation of Refs. 32) and 35) is that 
we improve the calculation of the angle part of the phase-space trace integral for 
the orbits with smaller degeneracies. Taking the elliptic billiard as an example, we 
have applied the ISPM to the integration over the angle variable for short diametric 
orbits, and derived the improved trace formula which is continuous through all bi- 
furcation points including the circular limit and the separatrix. We have then shown 
that significant enhancements of the shell effect in level densities and shell structure 
energies occur at deformations near the bifurcation points. Away from the bifurca- 
tion points, our result reduces to the extended Gutzwiller trace formula,^)' ^^"^'^^ and 
for the leading-order families of periodic orbits, it is identical to that of Berry and 
Tabor.'"^) 

The major purpose of this paper is to extend our semiclassical ISPM to the three- 
dimensional (3D) spheroidal cavity, which may be taken as a simple (highly ideal- 
ized) model for a heavy deformed nucleus^)' -"^^^ or a deformed metallic cluster, -^^^ 
and to specify the role of periodic orbit bifurcations in the shell structure responsi- 
ble for superdeformations. Although the spheroidal cavity is integrable, it exhibits 
all the difficulties mentioned above (i.e., symmetry breaking and bifurcations) and 
therefore gives rise to an exemplary case study of a non-trivial 3D system. We apply 
the ISPM for the bifurcating orbits and succeed in reproducing the supcrdeformed 
shell structure by the POT, hereby observing a considerable enhancement of the 
shell effect near the bifurcation points. 



*' In this paper SSPM denotes the standard stationary phase method and its extension to 
continuous symmetries. ^^"^^'''^ 
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§2. Classical mechanics for the spheroidal cavity 

The semiclassical trace formulas for the oscihating part of the level density for 
the spheroidal cavity are determined by the characteristic properties of the classical 
periodic families.^-*' ^^'^^^^^^^ This section is an outlook of the definitions and solutions 
of the classical mechanics for the spheroidal cavity in line of Refs. 8), 9), 23) and 27). 
They will be used for the semiclassical derivations of the trace formulas improved at 
the bifurcation points. We shall pay special attention to the 3D periodic orbits which 
emerge through bifurcations and play important roles as the semiclassical origin of 
super deformed shell structure.^^'^^^'^^^ 

2.1. General periodic- orbit formalism 

We characterize the spheroid by its ratio of semi-axes r] = b/a keeping its volume 
fixed, and consider the prolate case with ?? > 1, where the major axis coincides with 
the symmetry axis. We first transform the Cartesian coordinates (x, y, z) into the 
usual cylindrical coordinates {p,z,ip), where p = ya?+2/^, which are expressed in 
terms of the spheroidal coordinates (u, v, ip) 



p = C cosn sinhi>, z = C sin cosh t>, ( = \/W—a? (2-1) 

with 

-|<M<|, 0<'y<oo, 0<(^<27r. (2-2) 

The values of ib(^ define the positions of the foci of the spheroid lying on the z-axis. 
Taking into account the volume conservation condition a^b = E?, one has b = Rrj^/^ 
and a = Rri~^/^. As is well known, the Hamilton- Jacobi equations separate in the 
coordinates (u, v, ip) for the spheroidal cavity. 

In the Hamilton- Jacobi formalism, classical dynamics is determined by the par- 
tial actions. In the spheroidal coordinates these are given by 



lu = ^ [ du^l ai — sin^u ^ — , (2-3a) 



IT .l^uc V "^OS^ ^ 



= ^ / ' dvjcosh^v-ai (2-3b) 

TT V sinh V 

= \lz\ = pCV^, (2-3c) 

where Iz is the projection of the angular momentum onto the symmetry axis, and 
p = V^me, m is the particle mass. In Eqs. (2-3) we introduced new "action" variables 
o"! and (T2 related to the turning points —Uc, Uc and Vc, Vb along the trajectory in 
the {u,v) coordinates; u = Uc and v = Vc are the (hyperbolic and elliptic) caustic 
surfaces, 



coshvc = ^ 2^''" ~^ ~^ 



1 1 11/2^/2 



^(l-ai)2 + f72 



(2-4a) 
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-(l-(7l)2 + CJ2 



1/2 ■ 



1/2 



(2-4b) 



and V = V}j IS the spheroid boundary, given by coshv^ = r]/ \/rf — 1. The condition 
that the kinetic energy must be positive determine the hmits for the variables o"i 
and (72, 



172 < CTi < 



= < (72 < 



r/^ — 1 



0-2 (??^ - l) = (7 



+ 

1 ' 



= a A 



(2-5) 



These inequahties together with 27r intervals for the corresponding angle variables 
determine the tori of the classically accessible motion with the boundaries uf{a2) 
and (7^. 

According to Eqs. (2-3), the particle energy £ is a function of only the action 
variables and /<^, e = H {Iu,Iy,I^) due to integrability of the system un- 

der consideration. These relations define the partial frequencies w^, and 
through ujj = dH/dlj. The periodicity conditions for the classical trajectories are 
significantly simplified in terms of the partial frequencies cuj. Introducing the new 
variables k and 9, 



sm Uc 

cosh Vr 



arcsm 



in I 



/ coshi 



cosh V}) 



(2-6) 



along with the energy e instead of the partial actions /„, ly and (or ai and (72), 
they read 



1^ 



[(-(r)('-f1" 



(2- 7a) 



F(k) 



1 



1,2, 



F(k) 



(2-7b) 



Here, n^, n^, and are co-prime integers, n„ 
1,2, • • • , and k = ■sjrp' — l/{r] sin^). F and 11 are elliptic integrals of 1st and 3rd 
kinds (see Appendix A for their definitions). The periodicity condition (2-7) re- 
lates k((7i,(72) and 9{ai,a2) for a given periodic orbit /? to the integers ra„, n„ 
and n^, which together with the number of repetitions M define this orbit; i.e., 
f3 = M{ny,n^,nu). 



2.2. Three-dimensional periodic orbits 

The 3D periodic orbits (3DP0) M(nv,ni^,nu) form two-parameter (/C 



2) 



families for a given energy e since the number /C of free continuous parameters 
specifying an orbit with fixed energy and the same action is two.^^'^^ The condition 
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(a)ri=1.618... (b) ti=1.732... 




Fig. 1. The triangle of the classically accessible region determined by Eq. (2-5) is indicated by white 
lines in the (cti, (T2) plane at bifurcation deformations (a) 77 = 1.618 . . . and (b) r; — VS. The red 
and blue dots with circles indicate the 3DPO stationary points inside (really existing 3DPO) and 
outside ("ghost" 3DP0) of this triangle region, respectively. Several examples of the stationary 
points are indicated: on the (T2 = side, the short 2DP0 (elliptic triangle, square, and hyperbolic 
"butterfly"); on the (J2 = cri side, the short EQPO (triangle, square, star, and diameter (black 
crossed circle)). The long diameter (separatrix) is located at (ai — 1,(72 ~ 0). The color and 
the contour curves indicate (in unit of pC^/n) the curvature Kn defined by Eq. (3-14). 



for 3DP0 is the existence of real roots {k, 
deformation rj = ?7bif given by 



of Eq. (2-7). They appear at the 



Vhii 



sin(7^ny,/n^,) 
sin(7rnu/n„) ' 



1, 2, • • • , fly ^ ^ri^ -\- 1, Tiu^ — 2, 3, ' 



(2- 



where k = and 9 = tt{1 — 2nu/ny)/2, and exist for larger deformation rj > T^^if- 
These roots determine the caustics (the spheroid v = Vc and the hyperboloids 
u = ±Uc) of the periodic orbit M{nv,nip,nu) through Eq. (2-6). These caustics 
are confocal to the boundary of the spheroid v = Vh. 

Figure 1 shows the stationary points corresponding to the 3DP0 for two bifurca- 
tion points T/bif given by (2-8). The physical tori region (2-5) in the variables ctj is the 
triangle. At r/bif = 1.618 . . . (Fig. la), the stationary point for the 3DP0 (5, 2, 1) co- 
incides with that for the star-shaped (5, 2) orbit in the equatorial plane (discussed be- 
low) lying on the boundary with (T2 = <7i , and moves toward inside of the physical tori 
region for larger deformations. At ?7bif = 
\/3 (Fig. lb), the stationary point for 
the 3DP0 (6,2,1) lies at the bound- 
ary side and coincides with that for tri- 
angular orbits in the equatorial plane. 
At these bifurcation deformations, the 
lengths of the 3DP0 (5, 2, 1) and (6, 2, 1) 
coincide with those of the star (5, 2) 
and the doubly repeated triangle 2(3, 1), 

Fig. 2. Short 3D periodic orbits (5,2,1) and 
(6,2,1) bifurcated from the equatorial plane 
orbits (5,2) and 2(3,1), respectively. Their 
projections on the equatorial plane are also 
represented by thick-dashed lines. 
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Fig. 3. Some short moridian-planc orbits in the prolate spheroidal cavity: From the left to the 
right; the isolated long diameter (2, 1, 1), the elliptic triangular (3, 1, 1), the elliptic rhomboidal 
(4, 1, 1), the hyperbolic "butterfly" (4, 2, 1). 



respectively. Figure 2 illustrates some 
short 3DP0, and their projections onto 
the equatorial plane which remind us of 
the parent equatorial orbits. 

2.3. Orbits in the meridian plane 

Equation (2-7) have partial solutions for k((Ti, (J2) and (J2) which correspond 
to the separate families of orbits, i.e. two-dimensional periodic orbits (2DP0), in the 
meridian planes (containing the symmetry axis z) and in the equatorial plane. First, 
we consider the special solutions of Eq. (2-7) corresponding to the two-parametric 
(/C = 2) 2DP0 families in the meridian plane.^^'^^ For these orbits, (72 = and ai is 
in the regions 

2 

< (Ji < 1, 1< (71 < (2-9) 

T]^ — 1 

for the hyperbolic 2DP0 (with hyperbolic caustics u = iuc) and the elliptic 2DP0 
(with elliptic caustics v = Vc), respectively. The periodicity condition (2-7b) becomes 
the identity ujtp/cuu = 1 {n<f = 1, = 1), and is fixed by 

e = eh = arcsin [^^^^^ ' & = Oe = arcsin > (2-10) 

for the hyperbolic and elliptic 2DP0, respectively. For k, we only have the condition 
Eq. (2- 7a). This k determines ai and thus /„ and 1^ {lip = since (72 = 0) through 

^ 1 , , 

K = Kh = VOT) K = Ke = —=, (2-11) 



for the hyperbolic and elliptic orbits, respectively. 

Some examples of the hyperbolic and elliptic orbits lying along the triangular 
boundary side C72 = arc indicated in Fig. 1; see also their geometrical illustrations 
in Fig. 3. The hyperbolic and elliptic tori parts are separated by the separatrix point 
((7i = 1,(72 = 0) related to the long diameter (see below). Another end point of the 
hyperbolic tori coincides with the stationary point {ai = (72 = 0) for the diametric 



Symmetry Breaking and Bifurcations in the Periodic Orbit Theory: II 



7 




orbit in the equatorial plane. We can think of these hyperbolic and elliptic orbits 
as being periodic in the plane f = const and we call them "meridian-plane periodic 
orbits" . 

For the elliptic case, the solution k of Eq. (2-7) with 9 = 6e{n) exists for any 
Uu = 1,2, ■ ■ ■ and n„ > 2nu + 1, {n^p = n„) at any deformation rj > 1. Examples are 
the triangles ( u — 1), the rhomboids (4, 1, 1) and the star-shaped 

orbits (5, 2, 2) as one-parameter families in the meridian plane. The root k found 
from Eq. (2-7) gives the elliptic caustics with Uc = 7r/2 in Eq. (2-6) and the semi-axes 
ttc = CVl — k'^/k and be = C,/k,- 

For the hyperbolic case, the solutions k can be found for n„ = 1, 2, 3, • • • and 
even {n^ > 2{nu + 1)). In Fig. 1 the "butterfly" orbit (4,2,1) is indicated as an 
example. The families of these orbits appear for rj > %if with 



^if 



sm 



7rn„ 



nr. 



(2-12) 



This is the deformation where the diametric orbits M(2, 1) with M > 2 in the equa- 
torial plane bifurcate and the hyperbolic orbits emerge from them. Their hyperbolic 
caustics are expressed in terms of the root k of Eq. (2-7) and Eq. (2-6) with Vc = 0. 
The parameters Uc and be of these caustics are given by Uc = CVl — and be = (i^- 

2.4. Orbits in the equatorial plane 

In the equatorial plane with z = 0, the separate families of regular polygons 
and diameters are the same as in a circular billiard"^) of radius a. The restriction 
z = decreases the number /C to /C = 1. This single parameter corresponds to 
the angle of rotation of the polygons and the diameters about the symmetry axis z. 
Figure 4 illustrates the most important (shortest) equatorial-plane periodic orbits 
(EQPO); the diameters M{n,u = 2,n^ = 1), triangles M(3, 1), squares M(4, 1) and 
star-shaped orbits M(5, 2). They satisfy, from inequalities (2-5), 

<Jl =02, < (72 < (2-13) 

r/^ — 1 

Therefore their stationary points lie along the (72 = <ti side in the triangle, as indi- 
cated in Fig. 1. 

The caustic parameters Uc and Vc for these families are defined by Uc = and 
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Vc = arcsinh[acos(7rn(p/n,,)/C]- The solutions of Eq. (2-7) for these orbits are k = 



2.5. Diametric orbits along the symmetry axis 

In the spheroidal cavity, there is also a diametric orbit along the 2:-axis (see 
Fig. 3). It is isolated (/C = 0), since we have two additional restrictions x = and 
y = decreasing /C by one unit more than in the previous case. The solution of 
Eq. (2-7) for this orbit is k = 1 and 6 = arcsin(y^?72 — l/rj). Its stationary point 
coincides with the separatrix values (di = 1,(T2 = 0) corresponding formally to the 
caustic parameters {uc = 7r/2, Vc = 0); see the circle point with the vertical diameter 
in Fig. 1. (In Fig. lb, it is very close to the stationary point for the elliptic orbits 
(3,1,1) in the meridian plane, which lies slightly on the right along the (72 = side.) 

2.6. Bifurcations 

At the deformations r/bif given by Eq. (2-8), the EQPO M{nv,nip) bifurcate 
and the 3DP0 or the hyperbolic 2DP0 M{nv,nip,nu) emerge. We encounter the 
breaking-of-symmetry problem at these bifurcation points since the degeneracy (sym- 
metry) parameter /C changes there, for instance, from /C = 1 of the EQPO to /C = 2 
of the 3DP0 through them. Before the bifurcations (77 < rybif)) the stationary points 
(Tj of the 3DP0 and the hyperbolic 2DP0 are situated outside of the triangular tori 
region (2-5), and give rise to complex (k, 9) and complex caustics. Such formal orbits 
are called "complex" or "ghost" orbits. They cross the ai = U2 boundary through 
the stationary point of the EQPO at bifurcations [rj = r/bif) and then moves into 
the triangular tori region for larger 77. In Fig. 1 are also indicated the stationary 
points for the 3DP0 lying outside the physical tori region ((6,2,1) in Fig. la, (7,2,1) 
and (8,2,1) in Fig. lb). The equatorial diameters M(2, 1) correspond to the limiting 
case, o"! = CJ2 = 0. They bifurcate into themselves (/C = 1) and the hyperbolic 2DP0 
(2M, M, 1) in the meridian plane (/C = 2) at the deformations given by Eq. (2-12). 

The spherical limit is a special bifurcation point. Namely, the planar regular 
polygons (/C = 3) and diameters (/C = 2) bifurcate into the meridian 2DP0 (/C = 2), 
EQPO (/C = 1) and the isolated long diameter (/C = 0) for deformations r/ > 1. 
Another special point is the separatrix (ui = 1, C72 = 0) related to the long diameter. 
Near this point the complicated 3DP0, elliptic and hyperbolic 2DP0 having large 
values of {nu,n^) and nu/n^ close to 1/2 appear. Similar bifurcations of the 3DP0, 
EQPO and elliptic 2DP0 appear near other boundary values of cjj in the triangular 
tori on its "creeping" side ai = ai{a2) where some kinds of 3D "creeping" orbits 
having large values of riy but with finite and generally different n„ and appear. 
This is in analogy with the "creeping" singularities discussed for the elliptic orbits in 
the elliptic billiard^^) near the maximum value of ai, cj["^ = cosh^ = rj^ /{jf — 1), 
according to Eq. (A-6) at the right vertex in the "meridian-plane orbit" side 02 = 0. 
The 3DP0 with large number of the corners Uy and finite n„ = n^p approach the 
"creeping" elliptic orbits in the meridian plane. Another vertex corresponds to the 
creeping EQPO having large values of and n^p but for n^jn^p 1/2. 

The bifurcation point related to an appearance of "creeping" orbits cannot be 
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reached for any finite deformation. However, even for finite deformations like su- 
perdeformed shapes, the solutions for ai and (72 (related to the roots k and of 
the periodic-orbit conditions (2-7)) can be close to the "creeping" values of ai and 
CJ2 (related to their boundary values of (2-5)). In such cases we have to take into 
account such bifurcations in the trace formulas for the level density. The bifurcations 
of the 3D and 2D diameter- like orbits with Uv/uu close to 1/2 near the separatrix are 
rather long, however, so that they are not important for the shell effects discussed 
below. 

§3. Trace formulas for the prolate spheroid 

3.1. Phase-space trace formula in action-angle variables 

The level density g{e) is obtained from the Green function G{r,r';e) by taking 
the imaginary part of its trace: 

g{s) = Y,Sie-Sn) = ~lm J dr" j dr'G{r' ,r";e) 5{r" - r'), (3-1) 

where is the single-particle energy. Following Ref. 28), we apply now the 
Gutzwiller's trajectory expansion for the Green function G{r,r',e)}^''^^'^^^ After 
simple transformations,^^^ we obtain the phase-space trace formula in the action- 
angle variables {!,&), 

9sci{e) = ^ Re ^ I dO" J dl'S [e - H (/', 0')) 

IT 

(3-2) 



X exp 



1 (5„ (i',r, g + (r - J') • 0") - if 



where the sum is taken over all classical trajectories a, J = {/m,/^,,/,^} are the 
actions for the spheroidal cavity, = {Ou-,0v,&^} the conjugate angles, and Ua 
phases related with the Maslov indices. ^^^'^^^"^^^ The phase-space trace formula (3-2) 
is especially useful for intcgrable systems since the Hamiltonian H does not depend 
on the angle variables in this case, i.e., H = II{I). The action 

Sa{l',l",ta)=- j^^ dl-0{l) (3-3) 

is related to the standard definition 

1-0" 

Sa{0',@",e) = / dG 1(0) (3-4) 
J0' 

by the Legendre transformation 

Sa{@', 0", S) + I" ■ (0' - 0") = Sail', I", ta) + 0" • (I" - I'), (3-5) 

ta being the time for classical motion along the trajectory a. The phase Ua will be 
specified below. 
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3.2. Stationary phase method and classical degeneracy 

It should be emphasized that even for integrable systems the trace integral (3-2) 
is more general than the Poisson-sum trace formula which is the starting point of 
Refs. 5), 32), 35) for the scmiclassical derivations of the level density. These two trace 
formulas become identical when the phase of the exponent does not depend on the 
angle variables 0. In this case, the integral over angles in (3-2) gives simply (27r)" 
where n is the dimension of the system (77, = 3 for the spheroidal cavity), and the 
stationary condition for all angle variables are identities in the 2tt interval. This is 
true for the most degenerate classical orbits like the elliptic and hyperbolic 2DP0 
in the meridian plane and the 3DP0 with IC = n — 1 = 2. On the other hand, 
for orbits with smaller degeneracy like the EQPO (/C = 1) and the isolated long 
diameter (/C = 0), the exponent phase strongly depends on angles and possesses a 
definite stationary point. Therefore, we have to integrate over such angles by the 
ISPM in the same way as for the bifurcations of the isolated diameters in the elliptic 
billiard.28) 

3.3. Stationary phase conditions 

Due to the energy conserving (5-function, we can exactly take the integral over 
I'y in Eq. (3-2) and result in 

^^^'^'^ = (2^^^E / d&'uj do'^l de; I di'^ I 

(3-6) 



X exp 



'-{S^{l',l",t^) + {l"-l').@")-i^u^ 



The integration limits for 1^ and I^p are determined by their relations to the variables 
(ci, (T2) and by the boundaries given by Eq. (2-5). One of the trajectory ao in the 
sum (3-6) is the special one which corresponds to the smooth level density ^tf of the 
Thomas-Fermi model.^'^) For all other trajectories, we first write down the stationary 
phase condition for the action variables and I'^ : 

I^^MI^£M^ * -0';^0'^-0'; = 2^M„, (3-7a) 



dSaiI',I",ta)\ Q,,^ 27rM^, (3-7b) 



where M = {M.,^, M,^, M^) = M{nu,n^,n^) and Af arc integer numbers which indi- 
cate the numbers of revolution along the primitive periodic orbit /3. The superscript 
asterisk means that we take the quantities at the stationary point = I* and 
l!p = I^- We next use the Legendre transformation (3-5). Then, the stationary 
phase conditions with respect to angles {©u,&v,&(p) are given by 

d&' + 80' J =^ 
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In the following derivations we have to decide whether the stationary phase 
conditions (all or partially) given by Eqs. (3-7) and (3-8) are identities or equations 
for the specific stationary points. For this purpose we have to calculate separately the 
contributions from the most degenerate 3DP0, the 2DP0 families in the meridian 
plane (/C = 2) and those from orbits with smaller degeneracy like EQPO (/C = 1) and 
the isolated long diameter (/C = 0). The latter two kinds of orbit are different from 
the former two kinds in the above mentioned decision concerning the integration 
over angles 0. 

3.4. Three-dimensional orbits and meridian-plane orbits 

The most degenerate 3DP0 and the meridian-plane (elliptic and hyperbolic) 
2DP0 having the same values of action occupy some 3D finite areas between the 
corresponding caustic surfaces specified above. In this case the stationary phase con- 
ditions (3-8) for the integration over all angle variables 0„, 0y and 0^ are identities. 
The integrand does not depend on the angle variables and the result of the inte- 
gration is (27r)^. Since Eq. (3-8) is identically satisfied (the action does not depend 
on the angles like the Hamiltonian II{I)) we have the conservation of the action 
variables, I'^ = I'^ = lu and I'^ = I'^ = I^, along the classical trajectory a. The 
integrals over all O in Eq. (3-2) yield (27r)^ and we are left with the Poisson-sum 
trace formula,^-*' ^'^^ 



5sci(£) 



1 

¥ 



M 



Re^ / dIS{e-H{I)) exp 



M 



dl. 



exp 



— M -I -l-VM 

n 2 
27ri , _ _ . TT 

-—M-I-l-l^M 

n 2 



(3-9) 



It is convenient to transform the integration variables {Iu,I(p) to (cri,cr2) defined by 
Eq. (2-3), 



1 



5scU 



"2 da2 



, diu 

dO'l -r 



exp 



27ri , ^ ^ .IT 



(3-10) 



The integration limits are much simplified when written in terms of af^ (i = 1, 2) and 
form the triangular region shown in Figs. 1. We then integrate over crj expanding 
the exponent phase about the stationary point ai = a*, 

27r (M • /) = 5„ (/, t„) + (r - I) ■ 0" 

= Sp{e) + IY. - - ^j) + • • • > (3-11) 

i,3 



where S'/3(£) is the action along the periodic orbit /?, 

Sp{e) = 27rM [nJl + n„/„ (e, 7^, /;) + n^I*^] , 



(3-12) 
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and Iy{e, lu, Iif) the solution of the energy conservation equation e = H{Iu, ly, lip 



with respect to ly. Here the single prime index is omitted for simplicity. The J^- is 



the Jacobian stability factor with respect to crj along the energy surface, 

* i 

and K^j is the (2 x 2) curvature matrix of the energy surface taken at the stationary 
point (Tj = a* (at the periodic orbit /?), 

K^.= ^^-^^ I I ^^-^^ (3.14) 

daidaj Uy daidaj Uy daidaj 

See Appendix A for the explicit expressions of these curvatures. As we shall see 
below, the off-diagonal curvature K12 is non-zero for variables cjj. 

Then we use the ISPM, where we keep exact finite limits for integration over a,, 
and we finally obtain 

SgfLie) = - Re V exp (ikLp - i^up) , (3-15) 

12D/ ^0 ^ \ I J 

where eg = ^ /^rriE? {R'' = a?h due to the volume conservation condition). The 
sum runs over all two-parameter families of the 3DP0 or the meridian-plane (elliptic 
and hyperbolic) 2DP0, A^I^ is the amplitude for a 3DP0 or a 2DP0,*) 



4(2) _ 1 



■0 

dL, 



{2S} 47r {MnyR)^^a^\detKP\ [dai 



erf (Zf,Z+)erf [Z^^^t]- 

(3-16) 



Here, Lp represents "length" of the periodic orbit P 
i-KMuyp 



Lf3 = 

mujy 

= 2Mnybsme 



E{e, k) - E (k) + cot ^Vl-Ac2 sin2^ 



(3-17) 



where 9 and k are defined by the roots of periodic-orbit equations (2-7) {S^ = pLp for 
cavities). This "length" taken at the stationary points a* (the real positive roots 
of Eq. (2-7) through Eqs. (2-4) and (2-6)) inside the finite integration range (2-5) 
represents the true length of the corresponding periodic orbit (3. For other stationary 
points, the "length" is nothing else than the function (3-17) continued analytically 
outside the tori determined by (2-5). Wc introduced it formally instead of ujy by the 
equation uOy = 2TTpMny/mL(^ for convenience. In Eq. (3-16) wc introduced also the 

*^ The expression (3-16) is valid also for the 2DPO {02 = 0), since the product 172^^22 is finite 
for any (72 (see Appendix A). 
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generalized error function erf (Z ,Z~^) of the two complex arguments Z and Z~^, 

2 '•^"^ 



2 _ 2 

eif (^z~ , z~^) = / dze~^ = erf(2;"'") — erf 



(3-18) 



with erf (z) being the simple error function. The arguments of these error functions 
are given by 



± = ^J-iT:Mn,K(jh{af{al) - al) , 



Z^^ = J-iTrMny{detKf^/K^-^)/h {af - (Ta) , 



(3-19a) 
(3-19b) 



in terms of the finite limits af^ given by (2-5), and taken at the stationary point 
0-2 = 0-2. We note that, for the 3DP0 M{3t, t, 1) with t = 2,3, the curvature K^-^ 
is zero at any deformation. For such orbits, one should use 



2', 



/3± 



Z^^ = ^J-iTrMuyK^^/f^ 



K 



<^2 -(^2 + —g- [(^1 (0^2) - (^1) 



K: 



22 



(3-20a) 
(3-20b) 



in place of (3-19). The latter limits (3-20) are derived by transforming the integration 
variable a2 to C72 — {Ki2/K22){cri — crl). 

Let us consider the stationary points a* far from the bifurcation points. This 
means that they are located far from the integration limits. Accordingly, one can 
transform the generalized error functions to the complex Fresnel functions with the 
real limits and then extend the upper limit to 00 and the lower one to —00. In this 
way we obtain asymptotically the Berry-Tabor result of the standard POT,^) which 
is identical to the extended Gutzwiller's result^^ for the most degenerate (3D and 
meridian-plane) orbit families. 



A 



(2) 



1 



vr (Mn„i?)2yV|rdetl^ 



din 

dai 



(3-21) 



The constant part of the phase z^^ in Eq. (3-15), which is independent of rj and 
e, can be found by making use of the above asymptotic expression and applying the 
Maslov-Fedoryuk theory.^^)'^-*^)"^^) This theory relates the Maslov index /x/j with the 
number of turning and caustic points for the orbit family (3. For the 3DP0, the total 
asymptotic phase up is given by 

z^3D = H3B - ^630 + 2(Mn„ - 1), //3D = M(3n„ + 2n„). (3-22) 

Here denotes the Maslov index, the numbers of caustic and turning points tra- 
versed by the orbit, ep represents the difference of the numbers of positive and 
negative eigenvalues of curvature K^.*^ For the hyperbolic and elliptic meridian 

*^ Since the dimension of is 2, is written as eg = sign(_ftrf ) + sign(_R'2 ); where K'^ is the 
i-th eigenvalue of K'^. It is also calculated by e/3 = sign(iffj^) + sign(det if'^/Kf]^) for K^-^ ^ 0, and 
€/3 = sign(ir4) + sign(det K^^ /K^^) for K^^ ^ 0. Here, sign(a;) = ±1 for a; ^ and for a; = 0. 
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2DP0, one obtains 

z^2DH = /^2DH - ^e2DH + 2 (Mn^ - 1) , /i2DH = 2M (n^ + n„) (3-23) 

and 

t'2DE = A«2DE - ^e2DE + 2 (Mn„ - 1) , = 3Mn^ (3-24) 

respectively. Note that the total phase includes the argument of the complex ampli- 
tude (3-16), and depends on both deformation and energy. 

Near the bifurcation deformations, the stationary points a* are close to the 
boundary of the finite area (2-5). In such cases the asymptotics of the error functions 
are not good approximations, and we have to carry out the integration over ai in the 
calculation of the error functions in Eq. (3-16) exactly within the finite limits. One 
should also note that the contributions from "ghost" periodic orbits are important 
near the bifurcation points. It makes the trace formula continuous as function of r] 
at all bifurcations. 

When the stationary phase points a* are close to other boundaries of the tori, one 
has also to take the integrals with the finite limits; for instance, near the triangular 
side o"! = (J^{a2) where we have the "creeping" points for the 3DP0 inside the tori 
(2-5) and the meridian elliptic 2DP0 near the end point {a\ = af, G2 = 0) with a 
large number of the vertices — oo. Another sample of such special bifurcation 
point is the separatrix (fJi = 1,(J2 = 0) where 3DP0 and hyperbolic 2DP0 have a 
finite limit n^/n^ ^ 1/2 for ra^, ^ oo and n„ ^ oo. In this case the curvature Ku 
becomes infinite and the amplitude (3-16) approaches zero. Thus, to improve the 
trace formula near the bifurcations, we have to evaluate the generalized error integral 
erf{Z^ ,Z^^) (or corresponding complex Frcsnel functions^^^ ) in Eq. (3-16) within 
the finite limits Z^^ given by Eqs. (3-19) or (3-20). 

For the spheroidal cavity we have another bifurcation at the spherical limit 
where the "azimuthal" Jacobian ^i^d J\2 (^'l^) {^"2 ^ vanishes.^) This is 
the reason for the divergence of the standard POT result (3-21) in the spherical 
limit. Our improved trace formula (3-16) is finite in the spherical limit, because 

the "azimuthal" generalized error function erf(.22 ~, -Z2 ^) is proportional to J 22 
in this limit and thus this "azimuthal" Jacobian is exactly canceled by that coming 
from the denominator of Eq. (3-16). Thus, as shown in Ref. 9), the elliptic 2DP0 
term (/C=2) in tlic level density approach the spherical Balian-Bloch result for the 
most degenerate planar orbits with larger degeneracy (/C = 3), 



" t>l,q>2t V / V 



X Sin 



, „ / ttA 37r , TT 
2A;i?gsin — q - (t - l)Tr - - 

' J 2 4 



(3-25) 



where t = Muu and q = Muv Note that Eq. (3-25) can be derived directly from 
the phase-space trace formula (3-2) or from the Poisson-sum trace formula, both 
rewritten in terms of the spherical action-angle variables. 
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3.5. Equatorial-plane orbits 

We cannot apply the Poisson-sum trace formula (3-9) for equatorial-plane orbits, 
because, although the stationary-phase conditions for 0^ and 0" in Eq. (3-8) are 
identities, it is not the case for the angle variable &'^. We thus apply the ISPM for 
the integration over 0'^. 

Coming back to Eq. (3-6) we transform the phase-space trace formula to new 
"parallel" {0"; and "perpendicular" {0", 0"; J^, variables as explained in Ap- 
pendix B for more general (integrable and non-integrable) systems. We then make 
the integration over the (I^, I'^p) variables in terms of the ISPM by transforming them 
into the cjj variables. Next, we consider the integration over the angle variable 0'^ 
by the ISPM, as there is the isolated stationary point 0* = or integer multiple of 
27r. We expand the exponent phase in a power series of 0'^ about 0* = 0, 

S^{l,l",t^) + {I"-I).0" 

= pLeq + J^^"^ - - ^1) + i^uf + • • • ' (3-26) 

where the stationary point cr* = (t| = u* is given by 

Leq is the length of the equatorial polygon with vertices and M rotations, and 
is given by 

Leq = 2Mn^i?sin^, ^ = 'Kntpjny. (3-28) 
In this way one finally obtains the amplitudes for EQPO, 

'^5eq(^) = r E 4q exp [i (klEQ - liyEo) } , (3-29) 
^0 EQ z yj 



see Appendix B for a detailed derivation. Here, Zf^ arc the limits given by Eqs. (3-19) 
or (3-20) for i = 1,2, and = 0,Z+ = 2:+ from Eq. (B-19). The latter is related 
to the finite limits < 0„ < 7r/2 for the angle 0u in the trace integration in 
Eq. (3-6), taking into account explicitly the factor 4 due to the time-reversal and 
spatial symmetries. 

For the total asymptotic phase z^eq, one finds 

z^Q = Meq + ^, /XEQ = 3Mn^, (3-31) 

where //eq is the Maslov index. We calculated this phase using the Maslov-Fedoryuk 
theory^^^ at a point asymptotically far from the bifurcations. Note that the total 



sin^ 
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Table I. Bifurcation points of some short periodic orbits. 



periodic orbit 


»?bif 


periodic orbit 


rihii 


(4,2,1) 




(6,3,1) 


2 


(5,2,1) 


1.618... 


(7,3,1) 


2.247... 


(6,2,1) 




(8,3,1) 


2.414... 


(7,2,1) 


1.802... 


(9,3,1) 


2.532... 


(8,2,1) 


1.848... 







phase is defined as the sum of the asymptotic phase z^eq and the argument of the 
amphtude ^eq, Eq. (3-30), so that it depend on kR and r] through the complex 
arguments of the product of the error functions. In the derivations of Eq. (3-30) we 
have taken into account the off-diagonal curvature as in the previous subsection, but 
much smaller corrections due to the mixed derivatives of the action Sa with respect 
to O'l^ and CTj are neglected, taking cTj = cr* in Eq. (3-26). 

The bifurcation points are associated with zeros of the stability factor Pf^ and 
given by 

sin 6 

^if = ^^^5^' ri = l,2,--- ,M. 3-32 
sm {n(p/M) 

The bifurcation points most important for the superdeformed shell structure are 
listed in Table I. 

When the stationary points are located inside the finite integration region far 
from the ends, we transform the error functions in Eq. (3-30) into the Presnel ones 

and extend their arguments to ±00, except for the case when the lower limit is 
exactly zero. According to the definitions of the limit, Eqs. (3T9) and (B-19), for 
Zf, we have asymptotically +00 {i = 1,2,3), Z^ = Z^ — >• 0, Z2 for 

diameters and Z2 —00 for the other EQPO. Finally, we arrive at the standard 
Balian-Bloch formula^^ for the amplitude ^eQ' 



sm 



where /eq = 1 for the diameters and 2 for the other EQPO (erf(Z^,Z^) —>■ /eq in 
this limit). 

As seen from Eq. (3-33), there is a divergence at the bifurcation points where 
Fz"^ 0. We emphasize that our ISPM trace formula (3-29) has no such divergences. 
Indeed, the stability factor F^^ responsible for this divergence is canceled by F^^ 
from the upper limit Z^ , Eq. (B-19), of the last error function in Eq. (3-30), Z^ (x 

Ff^ , and we obtain the finite result at the bifurcation point: 



-^V^9;^^-(.r,.^)e.(.,,..^). ,3.3, 



It is very important to note that there is a local enhancement of the amplitude 
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(3-34) by a factor of order \/kR *^ near the bifurcation point. This enhancement is 
associated with a change of the degeneracy parameter K. by one unit locally near the 
bifurcation point and results from exactly carrying out one integration more than 
in the SSPM case. In general, any change of the degeneracy parameter /C by AlC is 
accompanied by the amplitude enhancement by (kR)'^^^^ because of the AlC extra 
exact integrations. These enhancement mechanism of the amplitude obtained in 
the ISPM is quite general and is independent of the specific choice of the potential 
shapes. 

We mention that a more general trace formula which can be applied also to 
non-integrable but axially symmetric systems can be derived from the phase-space 
trace formula (see Appendix B). 

The contribution of the equatorial diameters in Eq. (3-29) for deformations far 
from bifurcation points reduces to the Balian-Bloch result for the spherical diameters 
(/C = 2), 

Sgilie) = E 2^ sin{4MkR). (3-35) 

M 

The amplitude for the planar polygons in the equatorial plane vanish in the spherical 
limit (see Appendix B). Note that the contributions of the planar polygons in the 
spherical cavity, Eq. (3-25), are obtained as the limit of Eq. (3T6), for elliptic 
orbits in the meridian plane. 

3.6. Long diametric orbits and separatrics 

As mentioned in §2, the curvatures K^j become infinite near the separatrix 
(fji = 1,(72 = 0), sec Appendix C. This separatrix corresponds to the isolated 
long diameters (/C = 0) along the symmetry axis. Thus, for the derivation of their 
contributions to the trace formula, the expansion up to the second order in action- 
angle variables considered above fails like for the turning and caustic points in the 
usual phase space coordinates. However, we can apply the Maslov-Fedoryuk the- 
Qj,y39),4i)-43) similar way as the calculation of the Maslov indices associated 
with the turning and caustic points but with the use of the action-angle variables in 
place of the usual phase space variables. This is similar to the derivation of the long 
diametric term in the elliptic billiard. 

Starting from the phase space trace formula (3-6) we note that the spheroidal 
separatrix problem differs from the one for the elliptic billiard^^^ by the integrals 
over the two azimuthal variables 0'^ and which are additional to the integrals 
over 0'^ and I'^. We expand the phase of exponent in Eq. (3-6) with respect to the 
action and angle 0'^ about the stationary points 1^ = and an arbitrary 0* (for 
instance, 0* = 0), and take into account the third order terms, in a similar way as 
for the variables 0" and I'^ (sec Appendix C). Note that we consider here small 
deviations from the long diameters, and 0* determines the azimuthal angle of the 
final point r" of this trajectory near the symmetry axis. 

*^ The parameter of our semiclassical expansion is in practice ^JkLp ^oc \/kIi^ . It is actually 
large for 3D orbits (L/3 ~ IQR) associated with superdeformed shell structures in nuclei. 
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After the procedure explained in Appendix C, we obtain 



1/4 



2eoi2 "rfkR 

M 1=1 -\/ Co -c; 



x[A\{-w]) + iGi{-w] 



3' ^3' 



x[M{-wf,Zl,Zl)+iGi{-wf,Zl,Z+)] (3-36) 

(see Appendix C for notations used here). 

For finite deformations and sufficiently large kR, i.e. for large pC, oc kRyrf^—1, 
near the separatrix ai l,a2 —>■ 0, the incomplete Airy functions in this equation 
can be approximated by the complete ones. Thus, Eq. (3-36) reduces to the standard 
Gutzwiller's result for the isolated diameters,^)' 

with the length Lld(M) = 4Mb = AMrf'I'^R and the stability factor for long 
diameters given by Eq. (C-20). 

For the calculation of the asymptotic phase z^ld we use this asymptotic expres- 
sion and calculate the Maslov indices ;Uld by the Maslov-Fedoryuk theory,^^) 

VLT) = /XLD + 2, /XLD = 4M. (3-38) 

The additional phases, dependent on deformation and energy, come from the argu- 
ments of the complex exponents and Airy functions of the complex amplitude. 

In the spherical limit, both the upper and lower limits of the incomplete Airy 
functions in Eq. (3-36) approach zero and the angle integration has the finite limit 
7r/2 (see Appendix C). With this, the other factors ensure that the amplitude for 
long diameters becomes zero; namely, the long diametric contribution to the level 
density vanishes in the spherical limit. 



§4. Level density, shell energy and averaging 



4.1. Total level density 

The total semiclassical level density improved at the bifurcations can be written 
as a sum over all periodic orbit families in the spheroidal cavity considered in the 
previous section, 

Sgsciis) = Sg^^{s)+6g^i{s)+Sg'il{s) + 6gi^i{s) = ^^^^^^(e), (4-1) 

P 

where the first two terms represent the contributions from the most degenerate (/C = 
2) families of periodic orbits, the 3DP0 and the meridian-plane 2DP0, given by 
Eq. (3-15), the third term the EQPO given by Eq. (3-29), and the fourth term the 
long diameters given by Eq. (3-36). 
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4.2. Semiclassical shell energy 

The shell energy 6E can be expressed in terms of the oscillating part Sg^^-^ (e) of 



the semiclassical level density (4-1) as^^'^^'-^^^ 

E (^)''^5sd (^f), N = £deg{e). (4-2) 



— -» / r). \ 

SE 

' P 

Here, tp denotes the period for a particle moving with the Fermi energy ep along the 
periodic orbit 

tf, = MTf,= '^, (4-3) 

Tj3 being the primitive period (M = 1), M the number of repetitions, and the 
frequency. The Fermi energy e-p is determined by the second equation of (4-2) where 
A*" is the particle number. 

In the derivation of Eq. (4-2) we used an expansion of the amplitudes Ap{£) about 
the Fermi energy e = e-p- Although ^/^(e) are oscillating functions of the energy e (or 
fei?), we can apply such an expansion, because Ajj are much more smooth compared 
to the oscillations coming from the exponent function of kL^. The latter oscillations 
are responsible for the shell structure, while the oscillations of Aj^ merely lead to 
slight modulations with much smaller frequencies. 

Thus, the trace formula for 5E differs from that for 5g only by a factor {h/tp)"^ = 
{h?kp /mLjj)'^ near the Fermi surface, i.e. longer orbits are additionally suppressed 

periodic orbits. 



by a factor The semiclassical shell energy is therefore determined by short 



4.3. Average level density 

For the purpose of presentation of the level density improved at the bifurcations 
wc need to consider only an average level density, thus also avoiding the convergence 
problems that usually arise when one is interested in a full semiclassical quantization. 

The average level density is obtained by folding the level density with a Gaussian 
of width F: 

gr{e) = ^ / de'g{e')e-^'^^' . (4-4) 

The choice of the Gaussian form of the averaging function is immaterial and guided 
only by mathematical simplicity. 

Applying now the averaging procedure defined above to the semiclassical level 
density (4-1), one obtains^^'^^ 

SgrMe) = E '^^i? (-) --'"^'^ = E ^^i? (-) (4-5) 

13 13 

The latter equation is written specifically for cavity problems in terms of the orbit 
length (in units of a typical length scale R) and a dimensionless parameter 7, 

r = , (4-6) 
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where 7 is a dimensionless quantity for averaging with respect to kR. Thus, the 
averaging yields an exponential decrease of the amplitudes with increasing Lp and 
7. In Rcf. 9), the 7 is chosen to be 0.6. In this case, all longer orbits are strongly 
damped and only the short periodic orbits contribute to the oscillating part of the 
level density. For the study of the bifurcation phenomena in the superdeformed 
region, we need a significantly smaller value of 7. 

Finally, we can say that the higher the degeneracy of an orbit, the larger the 
volume occupied by the orbit family in the phase space, and the shorter its length, 
the more important is its contribution to the average density of states. 

§5. Quantum Spheroidal Cavity 

5.1. Oscillating level density 

We calculated the quantum spectrum by the spherical wave decomposition 
method^'^) in which wave functions are decomposed into the spherical waves as 

i^m{r) = Y^Ciji{kr)Yim{n). (5-1) 
I 

Here, m denotes the magnetic quantum number, and Yl' means that I is summed 
over even(odd) numbers for positive(negative) parity states, ji and Y;^ are the 
usual spherical Bessel function and spherical harmonics, respectively. The expansion 
coefHcients C/'s are determined so that the wave function (5-1) satisfies the Dirichlet 
boundary condition 

^m{r = Rin)) = (5-2) 

or equivalently, 

J dQYil,{Q)^rn{r = R{n)) = 0, ^l. (5-3) 

By inserting (5-1) into (5-3), one obtains the matrix equation 

Y^Bw{k)Ci,=Q, Bu'ik) = f dnYi*^{n)ji,ikR{n))Yi>^{n). (5-4) 
I' 

Truncating the summation / by sufficiently large number Ic, one can obtain the energy 
eigenvalue = h'^k^/2m by searching the roots 

dctB{kn)=0. (5-5) 

Figure 5 shows the energy level diagram for the prolate spheroidal cavity as functions 
of axis ratio rj > 1. In Fig. 6, we plot shell structure energy 

TV 

6E{N,ri) = ^en{v)-E{N,n) (5-6) 

n=l 
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Fig. 5. Single-particle spectrum for the spheroidal cavity as a function of axis ratio rj. Solid and 
dashed curves represent the positive and negative parity levels, respectively. 



as functions of i] and particle number N. As well as the strong shell effect at 
the spherical shape (r/ = 1), one clearly see a remarkable shell structure at the 
superdeformed shape (77 ~ 2). 

Next, we calculated the coarse-grained level density by the usual Strutinsky 
smoothing procedure by taking wave number k as smoothing variable: 

9,{k) =^-j^dk'RfM (^^^) (5-7) 

As a smoothing function fM{x), we took a Gaussian with M-th order curvature 
corrections 

/Mix) = -^e-^^Lj^/J^l^')' (5-8) 

where L'^i{z) represents a Lagucrrc polynomial. Eq. (4-4) corresponds to the case 
of M = 0. In the following, we took the order of curvature corrections M = 6 and 
smoothing width 7 = 2.5 with which we can nicely satisfy the plateau condition.^^^ 
The coarse-graining is also performed by the same smoothing function but with 
smaller 7. We define the oscillating part of the level density by subtracting the 
smooth part as 



Sgj{k) = gj{k) - gj{k). 



(5-9) 



22 



A.G. Magner, K. Arita, S.N. Fedotkin, K. Matsuyanagi 




Fig. 6. Shell structure energy 3E as a function of rj and N^^'"^, where A'' is the neutron (proton) 
number taking the spin-degeneracy factor two into account. Energies are counted in units of 
h^/2mR^ (~ 30A-^/^MeV). 



The left-hand side of Fig. 7 shows 6g'y{k) with 7 = 0.3, as functions of r/ and kR. 
One will note that a remarkable shell structure emerges at r/ ~ 2, corresponding to 
the superdeformed shape. 

Let us consider the mechanism of this strong shell effect. If a single orbit makes 
a dominant contribution to the periodic-orbit sum 

^5'sci(e) = X] cos(A;L/3 - vri/^/2), apik) = Af^/eo, (5-10) 

the major oscillating pattern in 5g should be determined by the phase factor of the 
dominant term. In that case, the positions of valley curves for Sg in the {rj, kR) 
plane are given by 

kLp-TTi^f3/2 = {2n + l)TT, (n = 0, 1,2,---). (5-11) 

The right-hand side of Fig. 7 shows the stationary action curves (5-11) for several pe- 
riodic orbits. Green solid curves represent the triangular orbit in the meridian plane. 
Other longer meridian orbits also make the same behavior but with smaller distances. 



1 1.5 2 2.5 3 3.5 4 
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meridian (3,1) 
equatorial (5,2) 
equatorial (7,3) 



3D (5,2,1) 
3D (7,3,1) 



. 7. Oscillating part Sg of the single-particle level density as a function of rj and kR (left-hand side) and stationary-action curves for several 
periodic orbits (right-hand side). Clear correspondence between the enhancement of the shell effect and the periodic-orbit bifurcations is seen. 
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Red dashed curves represent the star-shaped orbit with five vertices in the equatorial 
plane. It causes bifurcates at 77 = 1.618 . . . and the 3D orbit (5,2,1) appear (red solid 
curves). The sequence (n,2,l) (n = 5, 6, 7, • • • ) make the same behavior shifting the 
bifurcation points a little bit to larger r/. Comparing with the plot of quantum 5g, 
one notices a clear correspondence between the superdeformed shell structure and 
the bifurcation of above star-shaped orbits. One will also note the correspondence 
between the bifurcations of the equatorial-plane orbits (n,3) (n = 7, 8, 9, • • • ) with 
the hyperdeformed shell structure emerging at ~ 2.5. The significant shell energy 
gain at the superdeformed shape obtained in Fig. 6 is considered as a result of this 
strong shell effect in the level density. 

5.2. Fourier analysis of level density 

Fourier analysis is a useful tool to investigate quantum-classical correspondence 

in the level density.'') Due to the simple form of action integral Sfj = hkLp, one can 
easily Fourier transform the semiclassical level density Qsciik) with respect to k. Let 
us define the Fourier transform F[L) by 

F{L) = J dke-'''^g{k). (5-12) 

In actual numerical calculations, we multiply the integrand by a Gaussian truncation 
function as 

Fa{L) = f dke-^^^^'^\-'''^g{k). (5-13) 

v27r J 

Inserting the semiclassical level density (5-10), the Fourier transform is expressed as 
FfiL) = FML) + lX:e— ^/^a, exp [-1 (^) 



2 



(5-14) 



This is a function which has peaks at the lengths of classical periodic orbits L = Lp. 
On the other hand, we can calculate F{L) by inserting the quantum mechanical level 
density g{k) = Yln^i^ ~ ^n) as 

Ff'^L) = X] e-'^ik-^)^-^^-^ . (5-15) 



It should present successive peaks at orbit lengths L = Lp. Thus we can extract 
information on classical periodic orbits from the quantum spectrum. In the left-hand 
side of Fig. 8, we plot the Fourier transform (5-15) as a function of L and r]. In the 
right-hand side of Fig. 8, lengths of classical periodic orbits Lp{r]) are shown. Red 
curves represent the orbits M{ny, 2, 1) (n^ = 4, 5, 6, • • • ). We found strong Fourier 
peaks at 77 2± 2 corresponding to the periodic orbits (5,2,1), (6,2,1) and (7,2,1) just 
after the bifurcation points. We also found Fourier peaks at 77 ~ 2.5 corresponding 
to the periodic orbits (7,3,1) and (8,3,1) etc. Thus, we can conclude that those 
periodic orbit bifurcations play essential roles in the emergence of superdeformed 
and hyperdeformed shell structures. 




K3 
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2 4 6 8 



Fig. 9. Shell structure energies plotted as functions of N^^^ . Solid curves represent the exact shell 
structure energies. Dashed and dotted curves represent those calculated by using the coarse- 
grained level density g-y{k) with the smoothing width 7 = 0.3 and 0.6, respectively. 



5.3. Coarse-grained shell structure energy 

In order to prove that the shell structure at the superdeformed shape is due to 
the bifurcated orbits, we calculated the 'coarse-grained' shell energy defined by 

e{k)g-y{k)dk- / e{k)g;y{k)dk, (5-16) 

where the smoothed Fermi wave number k-p in each term is determined so that they 
satisfy the particle number condition 

/^f(7) rkvii) 
g^{k)dk= I g^{k)dk = N. (5-17) 

By coarse-graining with width 7, a shell structure of resolution AkR = 7 is extracted. 
Classical orbits relevant for such a structure are those with lengths 

L < W = (5-18) 

7 

Taking 7 = 0.6, contributions from periodic orbits with L^IQR are smeared out. 
Around the superdeformed shape, bifurcated orbits have lengths L ~ lOi? and those 
contributions are significantly weakened by smoothing with 7 = 0.6, and the major 
oscillating pattern of 5E should disappear if those bifurcated orbits are responsible 
for the superdeformed shell effect. In Fig. 9, the coarse-grained shell energies (5-16) 
calculated for rj = 1,2, 3, with 7 = 0.3 and 0.6, are compared with the exact shell 
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Fig. 10. (a) Semiclassical amplitudes I^IsdI for the 3DPO (5,2,1) and \Aeq\ for the EQPO (5,2), 
calculated at kR = 25 by the ISPM, are plotted as functions of the deformation parameter rj 
by thick and thin solid curves, respectively. They are compared with the SSPM amplitudes 
(dashed curves), (b) The same as (a) but for the 3DP0 (6,2,1) and the EQPO 2(3,1). 




Fig. 11. (a) ISPM3 amplitudes for the 3DPO (5,2,1), calculated at kR = 25, are shown as function 
of deformation r; by thick-solid curves. For comparison, the ISPM amplitudes and the results of 
exact integration in the Poisson-sum trace formula arc pfottcd with thin-solid and thick-dotted 
curves, respectively, (b) The same as (a) but for the 3DPO (6,2,1). 



structure energies. In the upper panel, one observes that the spherical shell struc- 
ture survives after smoothing with 7 = 0.6, indicating that the major structure is 
determined by orbits whose lengths are sufficiently shorter than lOi?. On the other 
hand, in the middle panel, one notices that the major oscillating pattern of the su- 
perdeformed shell structure is considerably broken after smoothing with 7 = 0.6. 
The same argument is valid also for rj = 3. This strongly supports the significance 
of bifurcated orbits for the superdeformed and hyperdeformed shell structures. 

§6. Enhancement of semiclassical amplitudes near the bifurcation points 

In this section, we present some results of the semiclassical ISPM calculation, 
which clearly show enhancement phenomena of the semiclassical amplitudes I^IsdI 
and I^eqI near the bifurcation points. 

Figure 10a shows the modulus of the complex amplitude A^d (Eq. (3-16)) for 
the 3DP0 (5,2,1) and ^eq (Eq. (3-30)) for the EQPO (5,2) as functions of the 
deformation parameter rf. They are compared with those of the SSPM. The SSPM 
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amplitude for the EQPO (5, 2) is divergent at the bifurcation deformation rybif = 
1.618 . . . , while the ISPM amplitude is finite and continuous through this bifurcation 
point with a rather sharp maximum at this point. This is due to a local change 
of the symmetry parameter /C from 1 to 2 at the bifurcation, and the associated 
enhancement of the amplitude is of the order VkR. As seen from Fig. 10a, the 
ISPM amplitude for the (5, 2, 1) is continuous through the bifurcation point and 
exhibit a remarkable enhancement a little on the right of it. It approaches the 
SSPM amplitude given by Eq. (3-21) away from the bifurcation point. The ISPM 
enhancement for the 3DP0 is also of the order V kR, because of the same change 
of the degeneracy parameter /C from 1 to 2 as in the case of the bifurcating EQPO. 
The same is true for the 3DP0 (6,2,1) and the EQPO 2(3,1) as shown in Fig. 10b. 

In Fig. 11, we consider corrections from the 3rd-ordcr terms in the expansion of 
the action about the stationary point. Here we incorporate the 3rd-order terms in 
the (Ti variable (ISPM3) which are expected to be important for the 3DP0 (6,2,1) 
whose curvature Kn is identically zero (see Appendix D). We also show results 
of exact integration in the Poisson-sum trace formula (3T0) (marked POISSON). 
One sees that the results of the ISPM3 for the (5, 2, 1) and (6, 2, 1) orbits are in good 
agreement with those of the ISPM in the most important regions near the bifurcations 
and on the right-hand sides of them. It is gratifying to sec that the ISPM and the 
ISPM3 amplitudes I^IsdI for (5, 2, 1) and (6, 2, 1) are also in good agreement with the 
results of exact integration in the Poisson-sum trace formula. With the 3rd-order 
corrections, excessive ghost orbit contributions in the ISPM (bumps in the ISPM 
amplitudes in the left-hand side of the bifurcation point) are removed and better 
agreement with the result of exact integration is obtained. Except for that, the 
corrections due to the 3rd-order terms are rather small, and good convergence is 
achieved up to the second-order terms. 

The amplitudes \Ap\ are slightly oscillating functions of kR. Since the period of 
this oscillation is much larger than that of the shell energy oscillation, one can use 
the expansion about the Fermi energy ep (or k^R) in the derivations of both the 
semiclassical ISPM shell energy SEg^i and the oscillating level density Sggd (3T5). 
Figure 12 shows the semiclassical amplitudes A^j) for the 3DP0 (5,2,1) and A-^q for 
the EQPO (5, 2) as functions of kR at rj = 1.618 . . . (top panel) and r] = 2 (bottom 
panel). In this figure, the semiclassical amplitudes for the 3DP0 (6,2,1) and 
A-EQ for the EQPO 2(3,1) are also plotted as functions of kR at the bifurcation point 
T] = \/3 (middle panel). We see that for r/ = 2 the amplitudes |^3d| for the 3DP0 
(5, 2, 1) and (6, 2, 1) become much larger than the amplitude |^eq| for the EQPO. 

§7. Comparison between quantum and semiclassical calculations 

In this section we present results of calculation of level densities and shell en- 
ergies with the use of the quantum Strutinsky method and the semiclassical ISPM, 
and make comparisons between the quantum and semiclassical calculations. In the 
quantum calculations, the averaging parameter 7 = 0.3 is used. 

Figure 13 shows oscillating level densities 6g for relatively small deformations; 
QM and ISPM denote the Sg obtained by the quantum Strutinsky method and the 
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Fig. 12. (a) Semiclassical amplitudes lAw] for the 3DP0 (5,2,1) and |^eq| for the EQPO (5,2) 
are plotted by bold and thin solid curves, respectively, as functions of kR at the bifurcation 
point r? = 1.618 .... (b) The same as (a) but for the 3DPO (6,2,1) and the EQPO 2(3,1) at the 
bifurcation point rj = 1.732 .... (c) Semiclassical amplitudes |^3d| for (5,2,1), (6,2,1) and |Abq| 
for (5,2) are plotted by thin-solid, thick-solid and dotted curves, respectively, as functions of kR 
at 7? = 2.0. 

semiclassical ISPM, respectively. For r] = 1.2 we obtain a good convergence of the 
periodic orbit sum (4-1) by taking into account the short elliptic 2DP0 with < 12, 

Hu = I, the short EQPO with the nicixiinurn vertex nuinber Pmax — (j^v)ma.yi — 

14, 

and the maximum winding number tmax = Mn^ = 1 (M = l,n^ = 1). The 
ISPM result is in good agreement with the quantum result. For the bifurcation 
point r) = of the butterfly orbit (4,2,1) and rj = 1.5 slightly on the right of 
it, the convergence of the periodic-orbit sum is attained by taking into account the 
contributions from the bifurcating orbits, (4, 2, 1) and the twice-repeated diameter 
2(2, 1) with t 

max — 2, in addition to the 2DP0 and the EQPO considered in the 

7/ = 1.2 case. 

Figure 14 presents the oscillating level densities for the bifurcation deformations; 
■q = 1.618 . . . for the EQPO (5, 2), r] = Vs for the EQPO 2(3, 1), and r/ = 2 for the 
triply repeated equatorial diameters 3(2, 1). It is interesting to compare this figure 
with Fig. 15 where some results of simplified semiclassical calculations are shown. 
In the top panel of Fig. 15, the SSPM is used instead of the ISPM. We see that the 
SSPM is a good approximation for rj = 1.2. In the middle and bottom panels, only 
bifurcating orbits are taken into account in the periodic-orbit sum: only the 3DP0 
(5, 2, 1), the EQPO (5, 2) and the butterfiy (4, 2, 1) are accounted for in the middle 
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panel, while only the 3DP0 (5,2,1), (6,2,1), (7, 2, 1),(8, 2, 1) in the bottom panel. 
By comparing with the corresponding ISPM results shown in Fig. 15, we see that, 
for = 1.618 .. . and 2, the major patterns of the oscillation are determined by these 
short 3DP0. 

Figures 16 and 17 show the shell energies which respectively correspond to the 
oscillating level densities shown in Figs. 13 and 14. Again, we see good agreement 
between the results of the semiclassical ISPM and the quantum calculations. For 
1] = 1.2, a good convergence is attained by including only the shortest elliptic 2DP0 
and EQPO, in the same way as for the level density 6g, see Ref. 9). For r] = \/2 
and 1.5, properties of the ISPM shell energies are similar to those considered for 
the elliptic billiard in Ref. 28). Now, let us more closely examine the bifurcation 
effects in the superdeformed region by comparing Fig. 17 with Fig. 18. In the top 
panel of Fig. 18, we show the ISPM result for ry = 1.618... in which only the 
bifurcating 3DP0 (5, 2, 1), the short EQPO (5, 2) and the hyperbolic 2DP0 (4, 2, 1) 
arc taken into account. In the middle panel of this figure, wc show the ISPM shell 
energies at rj = 1.732 . . . , calculated by taking into account only the 3DP0 (5, 2, 1), 
the bifurcating 3DP0 (6,2,1) and the EQPO 2(3,1). These comparisons clearly 
indicate that a few dominant periodic orbits determine the property of quantum shell 
structure at those bifurcation deformations. The bottom panel in this figure shows 
the dominating contributions of only a few shortest 3DP0 at ?7 = 2.0. Evidently, the 
short 3DP0 (5,2,1), (6,2,1), (7,2,1) and (8,2,1) determine the major oscillating 
pattern of the shell energy. Thus, we can say that they are responsible for the 
formation of the shell structure at large deformations around the superdeformed 




Fig. 13. Oscillating level densities evaluated by the semiclassical ISPM, and by the quantum me- 
chanics axe shown by dotted and solid curves, respectively, as functions of kR for several defor- 
mations r]. 
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Fig. 14. The same as Fig. 13 but for larger deformations. 
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Fig. 15. Comparison of the oscillating level densities calculated by quantum mechanics (solid 
curves) and those obtained by some specific somiclassical calculations (dotted curves); (a) the 
top panel shows a comparison with the SSPM result for r] = 1.2, (b) the middle panel shows 
the ISPM result in which only the bifurcating 3DP0 (5, 2, 1), the EQPO (5, 2), and the 2DPO 
butterfly (4,2, 1) arc taken into account for the POT sum in Eq. (3-15) for 77 — 1.618 . . . , (c) 
the bottom panel shows the ISPM result in which only the four shortest 3DPO are taken into 
account for 77 = 2.0. 
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Fig. 16. Semiclassical ISPM and quantum shell energies (in unit of eo) are plotted by dotted and 
solid curves, respectively, as functions of N^^^ . 
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Fig. 17. The same as Fig. 16 but for larger deformations. 



shape. These results of calculation are in good agreement with those obtained in 
Ref. 23) from the analysis of the length spectra (Fourier transforms of the quantum 
level densities). 
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Fig. 18. Comparison of quantum shell energies (solid curves) with those obtained with specific 
semiclassical calculations (dotted curves): (a) the top panel shows the ISPM result for rj = 
1.618..., where only the bifurcating orbits (5,2,1), (5,2) and (4,2,1) are taken into account, 
(b) the middle panel shows for rj = 1.732... the contributions of only the three orbits; the 
3DPO (5, 2, 1) and (6, 2, 1), and the EQPO 2(3, 1), (c) the bottom panel shows for r? = 2.0 the 
contributions of only the four shortest 3DPO to the ISPM sum. 



§8. Conclusion 

We have obtained an analytical trace formula for the 3D spheroidal cavity model, 
which is continuous through all critical deformations where bifurcations of periodic 
orbits occur. We find an enhancement of the amplitudes \Ap\ at deformations rj ~ 
1.6-2.0 due to bifurcations of 3D orbits from the short 2D orbits in the equatorial 
plane. The reason of this enhancement is quite general and independent of the 
specific potential shapes. We believe that this is an important mechanism which 
contributes to the stability of superdeformed systems, also in the formation of the 
second minimum related to the isometric states in nuclear fission. Our semiclassical 
analysis may therefore lead to a deeper understanding of shell structure effects in 
superdeformed fermionic systems - not only in nuclei or metal clusters but also, 
e.g., in deformed semiconductor quantum dots whose conductance and magnetic 
susceptibilities are significantly modified by shell effects. 
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Appendix A 

— Curvatures — 



A.l. Three-dimensional orbits 
The action is written as 

S = 2ttM (uyly + Uulu + n^I,fi) , 
where /„, ly and I,p are the partial actions. In a dimensionless form, 



I -^i 

TT 



I -^i 

TT 



TT 



where 



dz 



1 - ^ 

dz 



z"- 



1 



- ^2 )(z2 - 4), 



(z2-zi)(z2-4), 



lip = TT-v/o^- 

z± are related to the ctj variables by 

z'^ + Z^ = ai + 1, z'^z'i = (71 — (72- 

In terms of the elliptic integrals, (A-3) can be expressed as 

Iu = — [{z'i - l)F(fc) - a2U{z''_,k) + zlE{k)] , 
z+ 



(A-1) 
(A-2) 

(A-3a) 

(A-3b) 
(A-3c) 

(A-4) 
(A-5a) 



iv = — {{zl- z'i) [F {(p,k) -U{(p,n,k)]- zlE{(p,k)} + Zb sin (p, (A-5b) 
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with 



n = 



1-z 



'zr - z. 



2 ' 



If = arcsm \ 



z^-z 



+ 

2 ' 



zi, = cosh V}j = 



(A-6) 



Here, we used the standard definitions of the elhptic integrals of the first and the 
third kind*) 



F(99,A;) 



\/l — A;^ sin^ x 



(A- 7a) 



7o 



1 — A;^ sin x dx, 



(A-7b) 



U{(p,n,k) = / 

Jo (1- 



dx 



{1 — n sin^ x) — k"^ sin^ x 



(A-7c) 



and omit argument if = tt/2 for complete elliptic integrals. The action (A-1) is 
written as 



S = 2pCM 

The curvatures Kij of the energy surface s = H{ai, cj2, e) are defined as 

J^ij = — J^ii = 1 1- 



and the frequency ratios in Eq. (A-9) are given by 



di^/dai 



dluJi^ diu/dai 



(A-8) 



(A-9) 



(A- 10) 



dly _^ U)u diu 
da2 LOv d(T2 



(A-11) 



We used here the properties of Jacobians for the transformations from the variables 
{Iu,I(p) to (<7i,(72). For the first derivatives of the actions (A-3) with respect to ai 
and (72, one obtains 



1 



dai 



1 



22, 



Fi^,k), 



(A-12a) 



*^ The definitions of elliptic integrals (A-7) are related with those in Ref. 46) as F(9, k) = F{9\a) 
and n(0,n, /t) = n(n, ^|a) (K = sina). 
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^ = -— U{z^_,k) , 1^ = Cf F(^, k) + Cn n{ip, n, k) , (A-12b) 

(7(72 2;+ C(72 



with 



Cf = 



4-1 



1 



2z+{zl - 1) ' 

Cn = - 



zi-z^_ 



zi-z^_ 



2Z+CJ2 2z+(4 - l)(z2 - 1) 
For the second derivatives of these actions, one obtains 



(A- 13) 



dul 



d'lu 



2^{pW'^)-^^^)] (||-.^g)-||F(.) 



(A- 14a) 



2zlkj 



+ 



F((^ 



+ (A-14b) 



5(72 5(72 

+Cn 



F(^,fc) + CF 



9n A; dk 

1 dip dF{ip, k) dk \ dCn 



^in 5(72 



+ 



dk 5(72 



+ 



5(72 



(A- 14c) 

n((/?, n, fe) 



f dll{ip,n,k) dip dll{ip,n,k) dn dll{ip,n,k) dk 



V dip 



da' 



+ 



dn 5(72 



+ 



dk 5(72 



and 



5(7i5(72 24A;f 

d^iy ^ 1 

5(71 5(72 44A;f 



F(fc) + 



1 + fc^ 5F(fc) 
A; 5A; 



■ ((7i + l-2.;g)lang , 1 + fe^ 5F(y, fc) 
F {(p, k) : — TTT^, h 



A^zlA^ki 



k 



dk 



Here, 



Ax = \/l — sin^x, ki = — k'^, 



arcsin 



£± 



and 



54 _ 1 
5a7 ~ 2 



1 ± 



(71-1 



V(<T1- 1)2 +4(72 



1 

2 



1± 



4 + - 2' 



(A-14d) 



(A-14e) 



(A-14f) 



(A-15) 



(A-16) 
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aCp zl - 24 - 1 



2 I ^2 \ o-v2 /'-v2 ^2 \2 



da2 44(1 - z2 )2(^. _ ^ 

aCn ^ (72(34 + z^)-24(4- z^)^ 

5(T2 ~ Az%CtI{zI - Zl) 

d(f _ (2z2-(CTi + l))tan^ 
~ 2zlzlkfAl ' 

dn 0"! — 1 



da2 (1-4)2(4-4)' 



9(72 4 ~ ^- 

Derivatives of elliptic integrals are given by 

dll{ip, n, k) 1 



with 



dip (1 — nsin^ Lp)A^p ' 



dU.iip,n,k) 1 . ^, 

= - [n2i (<y2, n, k) - n(v?, n, k)\ , 



= ^ [ni3(99, n, /c) - U[ip, n, k)\ , 
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A. 2. Meridian-plane orbits 

For the meridian-plane orbits for which 1^ = (c72 = 0), the actions and 
defined by Eq. (2-3) can be simphfied. In the dimensionless form (A-2) one obtains 
for the eUiptic orbits 



Eie. 



1 

J ' 

1 



E 



+ 



Y^T/^ - - 1) 



(A-28a) 
(A-28b) 



Here we used the identity^^^ 



E((/?, k) — k^ sin cos ip/ \Jl — k'^ siv? 



(A-29) 



In this subsection, we omit the suffix "1" for the variable ai for shortness. For the 
hyperbolic orbits, 



in = 2 [E(V^) -{I -a) F(V^)] , 

I, = (1 - a) [F(V^) - YiOh, V^)] + nOh, V^) 



-E(V^) + 



vV'n^ - 1 



(A-30a) 



(A-30b) 



Equations (A-28) and (A-30) may be regarded as parametric equations in terms of 
the parameter a for the energy surface of the meridian-plane orbits, e{Iu, Ivi ^ip = 0), 
for its elliptic and hyperbolic parts, respectively. 

The curvature Kn of the energy curve for the meridian-plane orbits can be 
obtained by differentiating Eqs. (A-28) and (A-30) implicitly through the parameter 
a. In this way one obtains Eq. (3-13) with the dimensionless derivatives for the 
elliptic orbits 



diu 
da 



1 



1 



(A-31a) 



(A-31b) 



dl^ 
da 



2v^ 



■s/a 



■s/a 



(A-31c) 



9^ 
da^ 



n 



1 1 

o"' s/a 



no, 



1 1 



+ 



<7' V^J - (1 - (''^W 



. (A-31d) 



40 



A.G. Magner, K. Arita, S.N. Fedotkin, K. Matsuyanagi 



For the hyperbolic orbits, 



din 

da 



— [n(a,V^)-F(V^)], 



div 1 



d^iv 1 



da^ 4(7 
Thus, for eUiptic orbits 

1 



[UiOh, a, V^) - U{a, V^)z + F(V^) - ¥{6^, A)] • 



i^ii 



4V^ 



and for hyperbohc orbits 



i^ii = - 



1 



4(7 [ F(k) 
A. 3. Equatorial-plane orbits 

For the equatorial limit a2 = ai = a one obtains from (A-4) 

^2 = 0, 4 = ^ + 1- 



We thus obtain in this limit (k 0) 



din 



'PEQ 



dai 2^f^Tl' dai 2^^^TT' 
dzl \a/{a + lY 



and 



5(71 1 + 1) J ■ 

a'j^g ^ 7r(l - 2(7) 
5(72 8((7 + 1)5/2' 



8^ 
dai 



[a 



+ 1)5/2 



— {{2a- 1)(^EQ + 77 sin (2(^eq) 



2V5q^ [zg(l-(7)-((7 + l)] 



(^EQ = arcsm 
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Substituting (A-36) and (A-37) into (A-9) one finally obtains the equatorial limit 



kf? = ^>P"-'>;'- + ^' (A.38a) 



In the same way, one obtains 

zl{2 - a) + a{a + I) 



k^? = -by^-->^-y-^^J ^ (A.38b) 
8zla{a + irJzl-ia + l) 



Kf^ = + , (A.38C) 

The determinant of the curvature matrix for EQPO becomes 

det k^^ = ^4 (A-39) 

32z^a{a + 1)2 ^ ^ 

which is negative for any orbit and for any deformation rj > I. It shows that 
bifurcations of EQPO's occur only through the zeros of stability factor F^^. 

Appendix B 

Derivation of trace formula for the equatorial-plane orbits 



We start with the phase-space trace formula^) '^^^'^-'^^'^^^ 

xexp!^^[Sa{p',p",ta) + {p"-p')-q"] -^f^^a}, (B-l) 

where the sum runs over all trajectories a, q = qa{t, q",p') determined by the fixed 
initial momentum p' and the final coordinate q", H {q,p) is the classical Hamil- 

tonian, the phase related to the Maslov index, number of caustics and turning 
points. ■^^^'^''"^"^^^ The Sa ip',p",ta) is the action in the mixed phase-space represen- 
tation, 

fP" 

Sa{p',p",ta) =- dp-q{p), (B-2) 

related to the standard definition of the action Sa {q', q" , e), 

n" 

Sa{q',q",s) = dq-p{q) (B-3) 

Jq' 

by the Legendre transformations (the integration by parts), 

Sa {P', P", ta) = Sa («', q" , s) + {p' - p") q", (B-4) 



42 



A.G. Magner, K. Arita, S.N. Fedotkin, K. Matsuyanagi 



ta being the time for motion of the particle along the trajectory a. The J {p'±,p'±) in 
Eq. (B-1) is the Jacobian for the transformation from p'^ to p'_^. Here, we introduced 
the local system of the phase-space coordinates q = {qpq±} and p = {p\\,p±} 
splitting the vectors into the parallel and perpendicular components with respect to 
the trajectory a. 

For the equatorial-plane periodic orbits (EQPO) one of the perpendicular com- 
ponents q± and p± can be taken along the symmetry axis z, say z and pz, keeping for 
other perpendicular components the same suffix, q± and pj_. After the transforma- 
tion to this local phase-space coordinate system and integration over the "parallel" 
momentum = p = s/Ome by using the (5-function in Eq. (B-1), one obtains for the 
contribution from the EQPO (/C = 1) 



1/2 



exp\l[Sa{p',p",ta) + {p"-p')-q"] 



vr 



(B-5) 



where q^^ = dH/dp\\ = p/m is the velocity. In the spheroidal action-angle variables, 
^11 = Oy, p\\ = ly, gy = LUy, q± = Oy, = v?, p± = I^, z = 0u, Pz = lu, and we have 



1/2 



exp\'-[S^{r,l",Q + {l"-l').0"]- 



7T 



(B-6) 



We now perform the integrations using the expansion of the action Sa about the 
stationary points 



s.{i',i",Q + {i"-i')-0" 



(B-7) 



Here we omit the corrections associated with mixed derivatives of type d^S/d&dl 
for simplicity. J± is the Jacobian corresponding to the second variation of the action 
Sa with respect to the angle variable 0u, 



J 



EQ 



d'^S^ . d'^S^ . d'^Sa 



+ 2- 



+ 



dec dO'^dO'^ deCy^Q 



d0'u 



dl' dl" 
2 u _|_ '-'^11 



(B-8) 

This quantity can be expressed in terms of the curvatures K^'^ and the Gutzwiller's 
stability factor F^^, 



= 4 sin^ 



5^ 
1 



dV dl" 

2 u _|_ 



dO'' del' 



30'' 



EQ 



Mn„ arccos (l — 27/ ^ sin^ 



(B-9) 
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as 



J 



EQ 



EQ 



EQ 



In these equations we used simple identical Jacobian transformations 



(B-10) 



dl'u 



de'j. 



u / I\ 



d{eii'^) _de'i de':,di'^_ 
dii di'^dii 



The curvature is the quantity Ku defined in (B-13), evaluated at stationary 

point (Ti = (72 = cr* given by Eq. (3-27), and so on. 

The integrand of (B-6) does not depend on the angles (©i,, 6><^) and we obtain 
simply (27r)^ for the integration over these angle variables. We transform integration 
variables (/„,/<^) into {0x^02) to obtain simple integration limits, and integrate over 
(ci, C2) by the ISPM. In this way we obtain 



d{IuJip 



J±\det JEQ| 



X erf (.Ej^ , 2j_ ; 2^ , 2^ ; ^2 > ^2 



(B-ll) 



where 



erf ,x^;y ,y^;z ,2:+) = J dx J dy J 



dze-'^'-y'-'', 



(B-12) 



Note that the integration limits for the internal integrals over y and z in 
er{ (x~ ,x'^;y~ ,y~^; z~ , z^) in general depend on the variable of the next integra- 
tions, = y^{x) and z^ = z^{x,y). Here we define curvatures in the {Iu,Iip) 
variables as 



= 2Tr MuyKu, Ju, = 



dil 



= 2Tr MnyKtp, 



wp 



dludl^ 



2TrMnyKu^. 



Using (B-10) and relations 



det J = Jii J22 — J 



'12 



d{Iu,Iif 



d{<Ti,a2) 



{JuJ(p Juu:> )) 



(B-13) 



(B-14) 



irpa sm ( 



OJv 



irp 



ma sin (f) ' 



(B-15) 
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one finally obtains 

<^5eq(^) = J- XI (^^-^EQ - ^I'^eq) , (B-16) 



^0 EQ 



^EO = -\l ^-^^ erf (Z7,Z|;Zr,Z+;^-,^+) , (B-17) 

where -Leq represents length of the EQPO. The "triple" error function in Eq. (B-17) 
can be separated into the product of three standard error functions, 

erf {Zl^Zl- Z-,Z+; Z-,Z+) ^ erf {Z-,Z+) erf {Z-,Z+) erf {Z-,Z+) 

(B.18) 

by taking the limits at the stationary points for all deformations, except for a small 
region near the spherical shape. In this way we obtain the simple results (3-30). The 
arguments of the error functions are given by (3-19) or (3-20) for Zf (i = 1, 2) and 



^ 2^ 2h 16 V Mn„A;asin(^(7*(a* + l)detKEQ' ^ ^ 



The spherical limit is easily obtained by using the spherical action-angle variables 
{©0,&r,©ip;Ie,Ir:I(p}- In these variables 



1 / sin^ d) 

= .Mn/R,iF. (Z2,Zt;Z^,Z};Z-,Z}) , (6-20) 

where = n„ for the equatorial-plane orbits with (n„,ri<^), the invariant stability 
factor Fe = F^^ given by (B-9), 



' jp^^ j / \ 



(B-21) 



The quantities K^^ and K^^ 



\ l^S / EQ 

are the curvatures of the energy surface e = H{Iq, I^p) in the spherical coordinate 
system. In that system the maximum value of I,^ is equal to the absolute value of the 
classical angular momentum Ig, = ±7^, 7^ being the maximum value of |7e|, and 
= 0. We note that for the diametric orbits the stationary points Ig and 7* are 
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exactly zero and there are also specific integration limits in Eq. (B-20). In this case 
the internal integral over I^ within a small region can be evaluated approximately 
as 2/51, and one obtains for the "triple" error function 



We used here the fact that, in the spherical limit ^ the integral over Z0 can 
be approximated by the upper limit Z^^ given by Eq. (B-21). We omitted also the 
strong oscillating value of / dz^e"^^ at the upper limit since it vanishes after any 
small averaging over kR and equals 1 in this approximation. Wc also accounted 
for the fact that l/{TTpR) for the diameters; see Eq. (B-22) {(f) = tt/2 

for the diameters). Finally, the stability factor Fz is canceled and one obtains the 
Balian-Bloch result (3-33) for the contribution of the diametric orbits in the spherical 
cavity.^) 

For all other EQPO one has the stationary points I* = Ig ^ and I^ is identical 
to its maximum value Ig in the spherical limit. This is the reason why there is no next 
order (l/VkR) corrections to the Balian-Bloch trace formula for the contribution of 
the planar orbits with > 3. The latter comes from the spherical limit of the 
elliptic orbits in the meridian plane (3-16), see Ref. 9). 

Appendix C 

Separatrix 



Like for the case of the turning points,^^^'^^^ we first expand the exponent 
phase in Eq. (3-6) with respect to I^ : 



Sa {!', I", ta) - {I" - I') ■ 0" = c{l + cjx + c|a;2 + c!|x^ + . . . 



Here 



= tI+tIz + \z\ (C-1) 



^=\ilu- I'u) , (C-2) 



4 = 1 K - {I' - I'T ■ = {&',0",e) , (C-3) 



(C-4) 



' " ' 2pC,MhK ^00, ai ^ 1, (C-5) 



2 Kdii 



46 



A.G. Magner, K. Arita, S.N. Fedotkin, K. Matsuyanagi 

II fd^SaY 2Ti^h^M ( dk?,^ 



6 \dIL 



2^2 



dL 



^ < 0, ail ^ 1, 



(C-6) 



where the superscript asterisk indicates the value at = I" = I*. The asymptotic 
behavior of the constants c| near the separatrix ai 1 is found from 



log[(l + sin6l)/(l - sin 6*)] 



(ai-l)log3(ai-l) ' 
6 — 9h{r]) formally, see (2-6), 

dk^ 21og[(l + sin6')/(l-sin6')] 



<Tl 1, 



diu (((7i-l)log2(ai-l)) 



2 ' 



(71^1. 



(C-7) 



(C-8) 



The rightmost part of Eq. (C-1) is obtained by a linear transformation with some 
constants a and P, 



X = az + P, a 



= (3cS) 



-1/3 



(C-9) 



4 = (co - ciC2/(3c3) + 2ci/(27ci)) " , r| = a [ci - 0^/(803)] " . (C-10) 

Near the stationary point for cti — > 1, one obtains c\ and rj' —>■ — u)" with the 
positive quantity 



w" 



(3C3)4/3 



M log [(1 + sin 6l)/(l - sin 9)] (ai - 1) 



fllog((7l - 1) 



2/3 



(C-11) 



Using expansion (C-1) in Eq. (3-6) and taking the integral over 6>^' exactly (i.e., 
putting 27r for this integral), we obtain 



\ 



Ai {-wW , , Z+) + i Gi (-w;!' , , Z^'^ 



si^i:/''«:4,^[Ai(-»")+'0'(-"") 

a I III ^ C2 



where 



„ll r(c>-) 
''2 u 



(C-12) 



(C-13) 
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Ai(—w, zi,Z2) and Gi(— zi,Z2) are the incomplete Airy and Gairy functions defined 

by 



cos 
sin 



-wz + 



(C-14) 



Ai{—w) and Gi(— tf) are the corresponding standard complete functions, = 
Iu{o'^i^\ a^^^ ) is the "creeping" elliptic 2DP0 value defined in section 2. In the 
second equation of (C-12), we used the fact that for any finite deformation rj and 
large kR near the separatrix (cri — 1) 



0, 



Mlog[(l + sin0)/(l - sin6')]K' 
2(ai-l)2log^(ai-l) , 



1/3 



oo. 



(C-15) 



Using an analogous expansion of the action Tq in (C-12) with respect to the angle 
to the third order and a linear transformation like (C-9) one arrives at 



2£oTr^Rh 



, 1 (wllw^) 



1/4 



Ai(^-ii;ll) +zGi (^-wll)] ^Ai (^-w^,Zl,zf^ +iGi(^-w^,Zl,zf^ 



xexp\^[S:{l',l",e)-{l'-l"y-0n 



+■ 



2i 
3 L 



(^11)3/2 + (^±)3/2- 



TT 

« 2 f ' 



(C-16) 



where 



w 



( cl 
V(3C3)V3 



> 0, 



(C-17) 



W-' 



zt = 



TT 



3c^ 



1/3 



(C-18) 



+ 2 



+ 



F^J^ is the stability factor for the long diameters, 



xy 



-4sinh^ [Marccosh (2r]'^ - l)] , 



(C-19) 



(C-20) 
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4 



1 

Oh 



+ 3- 



'de'Jde'j. 



+ 



def. 



8 

u,a ^'^u,a 



< 0. 



(C-21) 



Note that according to (C-19) the quantity approaches zero near the separatrix 
((Ti — > 1) hke in the caustic case. This is the reason why we apply the Maslov- 
Fedoryuk theory"'^)' ^^-'"^^^ for the transformation of the integral over angle 0" from 
(C-12) to (C-16). The remaining two integrals over the azimuthal variables (/^ and 
O'^) can be calculated in a similar way as explained in the text. 

Divergence of the curvature Kip, Eq. (B-13), for the long diameters {ai ^ 1 , 
(72 0) can be easy seen from the following expression valid for any polygon orbit 
having a vertex on the symmetry axis, 



- -^0^ 



2rf 



1 + rp' tan^ ijj 



1 



(C-22) 



where Lq denotes the length of the side having a vertex at the pole, pQ the cylindrical 
coordinate of another end of this side, and ■0 the angle between this side and the 
symmetry axis. For the long diameters, Lq — 2hM, po — and -0 — 0, so that 
oo . 



Appendix D 

Derivation of the third- order term 



D.l. Third- order curvatures 

For the curvature ' which appears in the third-order terms in the expansion 
of the action Sj h with respect to ai , one obtains 



K 



(3) ^ d^Iy _|_ riu d^iu 
daf riy daf ' 



where 



daf 



2z\ 



dai 



6z 



dai 



— 3^4 



dzl d^l, 

dzl 



dai dai 



dzl 



= [n(^, k\ k) - F((^, k)] 4 - ^ F(<^, A;) + ^ ^ 



(D-l) 



(D.2) 
(D-3) 



B,, 



dk = 



dzl 



[n{k\k)-Y{k)\ dk-^nk), 



zl dk'' 



1 dzl 



dzl 



k"^ dai k'^ dai dai ' 



(D.4) 
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with the derivatives 



dai 



k 



+ TTT 



_ dFjk) 

ddk dk 

1 d'^zl d^zl 



k'^ d(j\ dai _ 

2(72 



^+[U{k^k)-F{k)] 

dF{k) dk dzl d'^zl 
dk 2z^ dai dai 



hdzi 
z^ dai 



F{k), 



dcrj [{ai - 1)2 + 4(72]'/' ' 



(D-5) 



dB^ 
dai 



dk 



dU{ip,k'^,k) dF{ip,k) ^ r_ 1 dz'i 



dai 



dai 



zt dai 



■+ 



Here 



dz'i dFjcp, k) J_ 
dai dai 



dU{k'^,k) k'^dk 



dai 

d_^_±d^\d^ d^ 
'dai Al dai ) dai + daf 



F{ip,k) 



dk 



2^2 



kjz 



+ (E(/c)-F(A:)) 



(D.6) 



(D.7) 



^ = lMk\k)-Fik)], 



(D-S) 



dU{(p,k^,k) 
dai 



kl 



I U{^, k\k) + ^ (E(^, k) - F{^, k)) 



(D-9) 



5^ " sin(2(^) \ i^zl - zlf 
dz'^ ( dz^_ 



\d^ U'-4)-^(^6-^-)J 



+ 2 



dai \dai 



12 2 \ ^Zj^ / 2 2 \ 



2cos(2^)('|i) , 



(D-ll) 



dai 



^ sm if + sin[2(p)- — 
z^ dai 



(D-12) 
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D.2. Stationary phase method with third-order expansions 

After the expansion of the action in the Poisson-sum trace formula (3-10) up to 
the second order with respect to (J2 and up to the third order with respect to cri, 
one obtains 



Re^ 



L 



Mn^Ryfal 



^ exp (zfcL, - 



/ da2 / 

Jo Jx- 



K 



dx exp ikC,MnyK22 [(^2 — (^2 + — ^i) 



where 



ci ^ 0, C2 = kCMuy^^^, C3 = ^kCMnyK[^\ 

K22 3 



X = ai — a1, x^ = af — al- 



After transformation from c72 to the new Z2 variable, 
^2 = \f--i 



kCMnyK22 0-2 - (72 + T^lo"! - 



and a linear transformation from x to 

x = qiz + q2, with qi = (Scs)"^/^, 

one obtains Eq. (3-15) with the ISPM3 amplitude 

Lr- ' 



92 



C2 

3C3' 



(D-13) 

(D-14) 
(D-15) 

(D-16) 
(D-17) 



ikC 



diu 
dai 



^Mn,R(kQMn,kf^)'n ^Mn,R'K22a*, 
X erf (2:2", 2:+ ) [Ai (-r, Z-,z+) + i Gi (-r, Z-, z+)] 



cxp I -ir^/^ 



Here Z2 is defined by Eq. (3-20b) and 



T = (3c3)-^/=^ 



3C3 



Cl 



a:"^ - g2 



In the limit ci -H' 



r = 



(3C3)V3 



(/cCMn^)^/3(detK/K22) 



(D-18) 



(D-19) 



(D-20) 



For finite curvatures far from the bifurcations, one can extend the limits of the 
Airy and Gairy functions {z- — ^ and 2;+ — ^ 00) and obtains the complete Airy 
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Ai(— r) and Gairy Gi(— r) functions. Then, using the asymptotics of these functions 
for large r oc {kKf/^ (large kR) 

[^\{-r) ~ ^1 ( -r3/2 + -1 (D.21) 

and of the erf-function in Eq. (D-18) one obtains the SSPM limit (3-21). 



